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ABSTRACT 


A model is developed which relates the structure of a 
combat unit to its expected effectiveness during a time inter- 
wal. Sensitivity to varvations Of structure is incorporated 
into the model by adapting effectiveness concepts from reli- 
۱ی‎ ۲٠۰ theory, and by constructing the attrition process as 
a continuous parameter Markov chain. The model is highly 
sensitive to the judgment of the user, due to parametric 
flexibility both in the attrition process, and effectiveness 
measurement. Two solution methods are given. Examples are 
employed to illustrate the model's use, and proposals for 


merinement are offered. 
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The purpose of any operational model is to aid a decision- 
Maker by relating his possible alternatives to their consequent 
results, thus increasing the knowledge upon which his decision 
is based. Therefore a model, in acdition to validly represent- 
ing the process or system involved, must be capable of dis- 
playing the contrast, in terms of meaningful criteria, of even 
subtle variations of choice. Sensitivity to variables outside 
the area of the current decision is not particularly valuable, 
and in fact may produce detrimental complication. In short, a 
model should ideally be a reflection of the problem of its user. 
It is for this reason that distinct decision problems often 
require the use of different models of a 001" system. 

This dependence of model upon problem is nowhere more 
apparent than in the study of combat dynamics. As an example 
a single battle situation may be modeled using game theory if 
the decision involves strategy, or using search theory if 
alternative fire distributions are being studied, or using 
Lanchester theory if weapon characteristics or tactical move- 
ments are being considered. Each modeling approach has been 
found useful for certain problem types but useless for others. 

The purpose of this thesis is to propose an alternate 
modeling viewpoint for those problems in which the structure 
Of a combat unit must be related to its effectiveness. The 


term structure is meant to denote the personnel and equipment 





of the unit, the roles of these elements, as well as the 
channels employed to receive and transmit information and 
instructions. As an example, a decision-maker who must 
choose for field test a small number of prototype tank killer 
teams from all the infinite organizational possibilities, 
faces a structural decision problem, and might profit from 
the model to be proposed. 

The importance of structure derives from the fact that 
when a leader or communicator becomes a casualty, his removal 
from the structure produces twin transformations in both the 
effectiveness and vulnerability of the unit. These transforma- 
tions are determined at least in part by the casualty's posi- 
Won in the organization. As an illustration consider a rifle 
Squad in defense whose squad leader has just become the first 
casualty. Eis function is assumed by a fire team leader, who 
must now receive and transmit orders, direct the fire of his 
squad and search his sector for targets. The effectiveness of 
the squad is reduced not only because one fewer weapon is 
firing but because a less highly trained leader is in charge. 
The vulnerability has increased for the team leader because he 
now must increase his exposure in order to direct the vhole 
Squad and observe its whole sector. The vulnerability of the 
whole squad is perhaps increased because it is being directed 
with less expertise. 

us q co. know the conc2tion Of a Unit xt rs not enough to 
know the number of casualties without knowing who they are. 


Perhaps the strongest testimony to this viewpoint is the 





universal military principle that enemy leaders and radio 
operators are the targets of highest priority in any type of 
engagement, and the simultaneous effort to conceal visual 
fans of their identity. If the validity of this approach 
has been justified the next step is to outline the model 
requirements. 

The decision-maker vill be assumed to have available to 
him a facility for military judgment which will include his 
own assessment of unit and individual vulnerability, and his 
own criteria for measuring unit effectiveness on the basis of 
its surviving structure. The mission of the model vill be to 


1) provide a framework in which the decision-maker's 


assessment of threat can be specifically quantified, 


2) translate the threat assessment into a stochastic 
description of the behavior of the candidate 


organization, 


3) provide a framework for expressing the decision- 
maker's effectiveness criteria as an algebraic 


PUNCELONe Of ties Unit econatt Lon, 


4) derive a stochastic description of effectiveness of 
the candidate organization, using the effectiveness 


2۶۹۰ 
Fortunately, the complete process of model construction is 
unnecessary, since there already exists an established field 
BEE Study which treats the concept of structure, The field is 
stochastic reliability theory, which has historically devel- 
oped in answer to the need for relating the design 7 complex 
Systems of mechanical and electronic components to system 


effectiveness. The comparison of a system oí components to 





a unit of fighting men may seem rash. However, through a 
process of adaptation, a model will be derived which may be 
considered more consistent while still taking advantage of 
Fellabi lity concepts. | 

Additional reliance on the theory of continuous parameter 
Markov processes will permit treatment of significant inter- 
dependencies in the attrition process. 

The resultant model is one with capacity for immediate 
application, as well as promising opportunities for improve- 


ment and refinement. 





7 eee Geo tie Ty CONCEPTS 


Reliability theory deals with collections of inter- 
dependent components. These collections are called systems. 

A component is any entity which can be described at any time 
aei ther functioning or failed, depending upon its ability 
to perform its assigned mission. Components are generally 
assumed to be in one of these two conditions. No partially 
failed conditions are considered, 

The system which is discussed in this thesis is an engaged 
combat unit whose components are included in one or more of 
EE following categories: 

1) personnel, 

2) weapons and equipment, 

3) rer Information zandseohkeundieatson. 

The number of components in the unit is represented by 
the variable n called the order of the unit. Each component 
mewrdentified with an integer in the set [1,...,n], called its 
index. 

The performance random variable of the ith component is 


defined as follows: 


X; Ey << ۵ If Ene 1 ات‎ component و‎ 6 
| failed condition at time t, 
= 1 if the ith component is in the 


functioning condition at time t. 
The component reliability of the ith component is defined 


as follows: 
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ri (t) = P[X; (t) = 1]. 


Since there are n components of interest, it is an obvious 


extension to define the performance random vector as follows: 


O E e). 


IE come Space ot X(t) eis the set of all binary n-vecters. 
Call this set 2 This set has N members where N = 2", 
It is often necessary to partition the outcome space se 


into m subsets, agressa where 1 < m < N, and define a 


me 
EcItion function by assigning’ a unique real number Di to 


all vectors in each subset aj, SO that 


E b: 7. ۶۶ vr X € ai. 


Each vector is a member of exactly one of the subsets. 
The first partition function of interest is the state 


Hu ction s(x), which is constructed by partitioning H. Into 


N subsets, each containing one vector. The value assigned to 
each one-member subset is one of the integers 1,...,N. 
Therefore 

s(x) - j if 7۶۲ X € à, and Di = j. 


Te is assumed here that in all cases 


Sl and 

 - N where 
, 4+4 707 and 
D NU UE 
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The composition of the state function with the perfor- 
Sg random vector results in definition of the state random 


variable 
S(t) = s(x(t)). 


The common practice in reliability theory is to consider 
the Nc system to be either in the functioning or failed 
condition at any time.  Predictably the system's condition is 
described by a binary random variable. However, when the sys- 
tem is a combat unit, such a measure is too coarse to be of 
value. Still, since this binary concept is of value in the 
adaptation process, it is introduced below. 

Mie structure function (x) of a systemwof ordesen 15 a 
partition function constructed by forming two subsets: 

1) a4, ت2۸‎ ٠7۰6 ۴  “-+ in Qr Such that if 


II Chen لا‎ SEEE  qnetionina, 


25 a», CHE ء٠‎ 1 ۶ D such that if 
Dai vico an the na legdweonaLleione 
Consistent with the use of binary indicators, 


(x) = 1(=b.) 0  " "+7 X € a, 


O (=b,) ir 57 X 6 a5. 
Composing the structure function with the performance 
mamaom vector results in definition of the structure random 


varlable 


Lr 





The system reliability is defined to be 
R(t) = P{d(t) = 1]. 


If more than one system is being considered then it is neces- 
sary to place a subscript on the structure function, struc- 
ture random variable and system reliability. 

A system is called coherent if its structure function 
satisfies three conditions: 

ول = لا 1 
L;‏ 


DD 

a 

ps 
I 


3) If x Ly, then g(x) < piy). 
Essentially a coherent system is one which fails if all its 
components fall, functions if all its components function, and 
Rose condition is not impaired by improving the condition of 
any of its components. 

DE conpeonenterscealled irrelevant irt its condition mever 
has any effect on the value of the structure function. 

Components are assumed to be of two varieties. Renewable 
components alternate between the functioning and failed condi- 
tions. For example, tactical radios are communication channels 
whose functioning is interrupted by atmospheric 8178611265701067 ۱ت‎ 
jamming.  Nonrenewable components never return to the function- 
Eug condition after entering the failed condition. Personnel 
are nonrenewable components of a unit. 

A common graphic method of illustrating the relationship 
of component performance to system ےا‎ is the use of 


the block diagram. Examples of these diagrams are shown in 


12 





Figures 1-4. The method is derived from the common practice 
Of representing electrical circuits. Each block corresponds 
to the like indexed component of the system. A potential 
flow is imagined, attempting to pass through the diagram 
along the connecting arcs. The flow is able to pass through 
٣٠٠٠٠٠۰67367. block only if the component is functioning. If 
the flow is able to traverse the diagram over any path then 
the system is functioning and its structure function equals 
one. Irrelevant components are represented by broken blocks 
totally unconnected to the system. 

Figure 1 shows the block diagram for the system of three 
components which functions if and only if all three components 
are functioning. This is called a series system. Figure 2 
shows a system which fails if and only if all its components 
p This is called a parallel system. Figure 3 shows a sys- 
tem which functions if and only if at least two of the compo- 
nents are functioning. This is called a two-out-of-three 
system. As in this case, a single component may be represented 
by more than one block, with the understanding that ali blocks 
with the same index are always in the same condition. Figure 
4 shows a system with components one and two in parallel, with 
component three irrelevant. It should be remembered that the 
block diagram shows only structural relationship, and implies 
nothing concerning actual location. 

A final concept, not unique to reliability, but useful 


here, is that of Laplace transformation. Given a function 
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Figure l. Series System 
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Figure 2. Parallel System 
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Figure 3.  Two-Out-of-Three System 
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Figure 4. Parallel System with Irrelevant Component 
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c(t), the Laplace transform of c(t) is defined as follows: 


o 


L {c(t)} = £&(s) = f exp(-st) c(t) dt, 
O 


E cunction of s. The function c(t) is called the inverse of 


the Laplace transform c(s). 
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TIT. MEASUREMENT OF EFFECTIVENESS 


In order to make any description of unit performance 
useful to a decision-maker, it is necessary to translate that 
description into some measure of effectiveness or goal attain- 
ment. Such a measure is useful only if it is consistent with 
the decision-maker's standards of value, and is in a form 
sufficiently simple to permit reasonable comparison of compet- 
ing alternatives. The method developed here, of constructing 
such a measure begins with the designation by the user of a 
relatively simple instantaneous indicator function based upon 
his intuitive value judgments. This function then serves as 
a basis for derivation of a one dimensional deterministic 
function, which measures unit quality over a variable time 
interval. | | 

The actual effectiveness of an engaged unit at a point in 
time is a very complex function of numerous variables, includ- 
ing friendly and enemy force levels, and many immeasurable 
physical and psychological factors. The simplifying assumption 
is made here that effectiveness depends only upon the vector of 
component performance and the length of time elapsed since the 
beginning of the engagement. The enemy performance is not 
considered directly, but the designation of an effectiveness 
function should reflect consideration of a hostile enemy envi- 
ronment, based on a previously determined scenario. The 
friendly unit is always assumed to have all components func- 


tioning at time zero. 
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The decision-maker's first step is to define a partition 
Bumetion on TE The set is partitioned into subsets ai, 
i-1,...,m, where vectors in the same subset are judged to 
represent equal effectiveness unit conditions.  Vectors in 
different subsets represent outcomes which are not equally 
effective. This partition is performed without reference to 
time. To each subset a; is assigned a dimensionless real 
number gi; called its absolute effectiveness level. The 
function thus defined is called the absolute effectiveness 


Function g(x) where 
g(x) = g; 2 ۶ ۱ 16 x ea 


and the partition is constructed so that 


Tt is assumed without serious limitation that 


7 0 و رہم 
g(1) = gv,‏ )2-2( 
Da SY then g(x) € g(y).‏ )3-3( 


These assumptions merely require that a completely destroyed 
unit be considered conpletely ineffective, that a completely 
undamaged unit be considered maximally effective, and that 
loss of a component never improve effectiveness. 

Since the vector members of 2, and the integers in the 
set {1,...,N} are placed in one-to-one correspondence by the 
deve function s(x), it is possible to designate a companion 
aunotion to g(xX) called the alternate effectiveness ۶۰77) 


definedton {N such thet 
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d(j) = 9, 


i i1 ۷ IE E ۳۱ E a; 


It follows that 
Slee) T 


This function is useful in later calculations. 

The construction of the actual effectiveness function is 
of extreme importance to the successful employment of the 
model. The user should take special care to assure that the 
absolute effectiveness levels truly represent his judgment, 
not only vith respect to ordinal relationships, but to ratio 
releticnships as well. 

A unit's real effectiveness depends upon how long it has 
been engaged. One factor to be considered is the change in 
enemy situation due to losses or reinforcement as provided in 
the scenario. Another consideration is that in missions of 
fixed length, such as delaying actions, casualties are less 
damaging, the later they occur. To allow for such trends, the 
۰٠۰٠٦ next step is to define a time variance function, f(t), 


and a resultant actual effectiveness function, 
h(x,t) - g(x) £(t). 


۱۳۱ a Specific application the variance effect 1s not con- 


Sidered, then 


f(t) 
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and therefore 


h(x, t) 
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۱۱ ۳۱ ۶۱ 0 OL G(x) can be expressed in another equiva- 
lent manner which shows the close relationship with reli- 
ability theory.  Esary, Proschan, and Walkup [Ref. 3] point 
Ew that a function such as gí(x), having a finite range can 


be written as 
34) g(x) = (9,-9)) GB) (x) + (95-95) Bo (x) + 
+ (9, 794.1) رر‎ (x) 
where the functions D, (x) are defined such that 


(3-5) Ø; (x) > 0 1 ۶ and Ol if g(x) < Sh 


Il 


(3-6) Ø; (x) 1 3270007 8 queue gi. 


It is now shown that each of the functions Ø. (x) is the struc- 
ture function of a unique binary coherent system composed of 
the components of the unit, in which at least one of the 
components is not irrelevant. Furthermore the systems are 
related in such a way that no system is functioning if any 
lower indexed system has failed, and no system is failed if 
any higher indexed system is functioning. 
Morem 1 

The functions Ø. (x) in (3-4) are structure functions of 


coherent systems of the n components of the parent system. 


Proof 
۳ mask be shown that for و ۰۰ و11‎ 
1) ¢,(0) = 0, 
2) 9, (1) > 


E ea) cy). 
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The proof is by contradiction. The above three points are 
each in turn assumed false, and a result, contradicting one 
of the assumptions made in defining g(x), is shown to follow, 


1) Assume that Ø. (0) z 0 for some i. Then it must be 

that g; (0) = 1, but then from (3-4) 

g(0) ES CE) < 0 

contradicting (3-1). “Thus Pp. (0) EIA o Sl 0. 
2) Assume that 2, (1) A 1 for Some +. Then it MUSE be 

chat Ø. (1) = 0, but then from (3-4) 

ES (9,797) + o... + (3. 784, 1) کت رس که‎ 
GONE E 

conbEadrctsuo (3-2).— Thus Ø. (1) = Ll, for iSl.. Sme 
۳۱۳ Assume that there exist vectors x and y in Gr such 

thgtox Sand that for some Ir Ø. (x) > Pi (Y). Since 

the function is binary, Ø. (X) = ], and Ø. (y) = 0. By 

wo), Ø. (X) = 1 implies that g(x) > gi. By Qe 

$. (y) = 0 implies that g(y) í g;. Therefore 

and X £ Y‏ رت 
۵۸ 8ہ 1 contradicting OE THOS ILE‏ 
$.GU < d;ly), for à = 1,..., ٤71‏ 


Theorem 2 


The functions Ø.: (x), i = l,..., Mel, are related as follows: 
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2) 


Proof 
The 


1) 


2) 


Thecrem 


For 


203 


1 
pa 
x 

۱ 
LA 


8 By (x) be hen D Lx) lysstor 3 


Il 
Il 


gef: Bj x) 0, then 9. (x) DATO: جا‎ e. 


Proof 15 agama by Contradiction. 


Let Ø. (x) LTO CE Vert ITE then the 
result is trivially true. If k > 1, assume that 
pe (x) = 0 for some j < k. By (3-6), Ø, (x) - ] implies 
E Ipe By (3-5), ۵ (x) = 0 implies that 
a) < gj: Since 9. * gy; Ve fol ove erat 

g(x) > g} < g(x), 
an obvious contradiction. Therefore A = l for 
ESI no kml. 
Let Ø, (x) `= on T e 
result holds trivially. If k < (m-1), assume that 
Pa (x) =) for Some | By 039509 Dy (x) = 0 implies 
that g(x) < g,. By (3-6), RES) = ] implies that 
Box > 3 Since he > ہو‎ it follows that 


g(x) < g4 * g(x), 


an Obvious contradiction. Therefore Es = 0, for 

TE +10 AN | 

3 

any vectors gx) 5 gi, ظ2‎ only if 203 سے‎ 0, and 


airs 
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Proof 


1) Let g(x) = رو‎ for some x. Then by (3-5), g(x) = 


gr 2.9. implies that 2; (x) ou DEV So XS s 
بت‎ implies that Bi )( = 1, 

2) Let g, (x) = 0, and g;_; (x) = l for some x. By (3-5) 
Bi (x) SS ns gi. By (3-6), 
g; q (x) < 11۲00001066 ۰) 91º Therefore it 


follows that 


but the only value in the range of the function which 
satisfies the inequality is g,. Therefore g(x) = g.. || 
Pete ultimate importance of this method of expressing g(x) 
is that the study of a very complicated multivalued absolute 
effectiveness function can be simplified to the study of a 
linear combination of simple dependent binary structure func- 
tions of the type commonly treated in classic reliability. 
It follows that the actual effectiveness function can be 
written as follows: 


MESE) = E(t) i{go-g,) 6) (x) + s t (S. ug. در‎ 


Thus far the functions introduced are deterministic with 
domain ° Since 2 is also the outcome space of the perfor- 
memes random vector, X(t), Sel of the functions can be 
composed with this vector to define analogous random variables. 
The absolute effectiveness random variable is defined by the 


relation 





no 


e - — Ce ” 


۱ > 2 
۲ 





or equivalently 


G(t) = (95794) 9, (t) TUO S I4 (d 794. 1) op ft). 
where 


9 (E) = 6, (X(t)), i = 1,..., m2, 


The actual effectiveness random variable is defined by the 
relation 
۲۱۳۱ ۲ ۱ (۵ 

These effectiveness random variables have the advantage 
that they are defined through a process which assures their 
consistency with the decision-maker's standards. Their pri- 
mary limitation is that they represent only instantaneous 
effectiveness. That is for every nonnegative value of t 
there is cefined a distinct random variable H(t). In Figure 
5, the bold horizontal line depicts a possible curve of out- 
comes of H(t), where f(t) is assumed to equal one for all non- 
negative t. 

Since combat engagements are extended over indefinite time 
p 5Nods, it is desirable to aggregate the effectiveness curve 
into a variable which measures accumulated effectiveness over 
an interval. To fill this need the value random variable is 
defined by the relation 


E € 
V(t) = J H(s) ds = ۱ G(s) f(s) ds. 


O O 
In Figure 5, the value random variable is represented by the 


Shisded area below the curve of H(t). 
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Example Outcome of Actual Effectiveness Random 


Variable and Value Random Variable 
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The value random variable is a measure of the type 
desired, but in general it will have a complicated distribu- 
Bon, too unwieldy to be of practical use in decision making. 
A relatively simple measure of effectiveness, which still 
preserves many advantages of the value random variable, is 
EU guuality function defined by the relation 

A ANA 
the mean or expected value of the value random variable. 

To compute the quality function in a specific case, note 

phat 


E 


E | H(u) du 


O 


q(t) 


te 
f em) du 


O 


It 


E 
f f (ú) E[G(u)] du, 


O 


where the interchange of expectation and integration in the 
Eng line is justified by application of Fubini's Theorem, 
since H(u) is a family of nonnegative random variables. The 
validity of this interchange is established by Royden [Ref. 
E]. 

There are two approaches for calculating q(t). The first 


merses from the fact that 


. 
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E[G(u)] = ) g(x) PIX(u) = x [X(0) = 1) 
A 
N 
= ) dali) piu), 
ci 
where 
A O OS 


is called the interval probability function from state one to 


ELLO 1, Thus the first result is that 


È N 
ene) | £ (u) | ) à) Pii LV E 
a 


Il 


O Em] È 
N t 

= ) dti) 1 f (u) Pj (2) du 
1=1 O 


Il 
ee 
O, 
E 
HH 
c 
e 


(3-7) 
1-1 
where 
t 
(3-8) I(i,t) - f E) pj; (uU) du. 
O 
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mic Second appreach follows from relation (3-4), 


mol 
E[G(u)] = EI 2 (9341795) *, Q1 
al 
mz 
z re oleae 
j=l 
meL 


m Se plo, (u) = 1 


Li. 
Il 
| 


T 
te 


= ee RCE R; (u) , 


C. 
Il 
pal 


where Bom) is the system reliability of the jth associated 
System, computed on the basis of all components functioning 
at time zero. And the second result is that 


16 m-1 


et) - | f(u) T gdo ۵ du 
O j=l 
11-1 t 
= SCI Re 
(9441 ga) S ; (u) du 

j=l O 

101-1 

31 

where 
E 
(3-10) 3ك‎ )3( = fe Bos Cu. 

O 


Now since the functions f(t) and g(x) or d(i) are assumed 


already defined by the decision-maker, all that is required for 
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calculation of the quality function is either set of 
functions pP) Ware... Pin (MS, or (R, (u),..., R _j(u)}. The 
determinetion of the functions is treated in (Section (VE 

The quality function is the primary measure of effec- 
tiveness proposed here. It is a one dimensional deterministic 
function, yet it is based upon consideration of the random- 
ness of unit performance. Since it is derived from functions 
defined by the decision-maker, it is consistent with his 
Wibiicit judgment. And since it is a function of time, it 
displays sensitivity to the length of engagement, 

In order to illustrate the use of the quality function in 
predicting unit effectiveness, this section will be closed with 
an example. 

Example l 


Consider a tank killer team consisting of the following 


Components: 
COMPONENT INDEX 
Team Leader 1 
Gunner 2 
Recoilless Rifle 3 ۲ 


Memes assumed that the unit is engaged in a large scale battle 
as part of a large force. The performance random vector for 
this unit, X(t), has an outcome space Vas consisting of the 
eight binary 3-vectors. Let the state function be defined as 


follows: 


2L 





5) )1,( , 1 1 ٭>ت-‎ ۶۵۹8 
SMA eee s{(0,1,0)} = 6 
s{(0,1,1)} = 3 Stio, 0, ۲ ۶ 7 
s((1,1,0)) = 4 s{(0,0,0)} = 8 


The decision-maker is now faced with the problem of 
Selo eructing the absolute effectiveness function g(x).  Sup- 
pose the hypothetical team is considered to have two missions. 
At any time the weapon is functioning, the primary mission is 
to detect enemy tanks and bring them under fire. If the 
weapon is in the failed condition, the secondary mission is to 
detect enemy tanks and direct the fire of adjacent teams 
against them. From experience, field test, or some other 
information source, the user determines that when both person- 
nel and the weapon are intact the team is capable of bringing 
effective fire on eight tanks. per hour. He decides to base 
his absolute effectiveness function on this fire capability, 
ENS first step is to set g{(1,1,1)} = 8. He feels that 
in most units the leader and gunner will have equal ability in 
firing the weapon. Therefore when either member becomes a 
casualty when the weapon is functioning, the decrease in effec- 
tiveness is about the same. The user estimates the reduced 
Capability to be about five tanks per hour, and defines 
g{(0,1,1)} = g{(1,0,1)} = 5. When the weapon fails, the team 
performs the secondary fire direction Mission. “The user dererz 
mines that both members together can detect and bring fire on 
as many as three tanks per hour. However, if the leader alone 


survives, he has a capability of two tanks per hour, while the 
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gunner alone has a capability of only one per hour because he 
has less experience and training than the leader in directing 
fire. In the two cases where both members are casualties the 
capability is, of course, zero. Thus the decision-maker 
Gefines his absolute effectiveness as follows: 

a prt te acl) eS gg = 8 

ie) go = 5 

g{(0,1,1)} = d(3) = g, = 5 

eta = d(4) = 23 ٦ 

Gel C0 TS g3 = 2 

SUD I ducc) s 1ہ‎ 

COICO OI 2 = g, = 0 

٥‏ ۹۷۹ ) :ئ00 
The only reguirement to be met in practice is that the‏ 
resulting function faithfully reflect the true values of the‏ 
decision-maker. Additionally, when several units are being‏ 
compared, a Separate function must be defined for each one,‏ 
and however they are constructed, it must be the case that‏ 
states of different units adjudged by the user to be equally‏ 
de c tive have the same value of g(x).‏ 

Femina cated IN (24), ORL GaN be expressed jac folte.<: 

g(x) = (1-0) g (x) + (2-1) go (x) + (3-2) $4,600 + (5-3) da lx) + 
(8-5) Ø. (xX) where the structure functions Ø. (x) are defined 


in the table below, 
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x B, (x) Ø, (x) 23 MES do (x) 
i51, 1) 1 1 Į 1 1 
B0. 1) 1 1 1 1 0 
80 1, 1( 1 il 1 1 0 
e, 0) 1 1 1 0 0 
B 0) 1 1 0 0 0 
i. 0) 1 0 0 0 0 
070,1) 0 0 0 0 0 
070, 0) 0 0 0 0 0 


The block diagrams of coherent systems associated with these 
Structure functions are shown in Figure 6. 

The decisicn-maker also wishes to define a time variance 
function for the tank killer team which enhances the actual 
effectiveness of every state the longer the EUM lasts. 
He feels that enemy losses and confusion make friendly fire- 
power more valuable later in the battle. He consequently 
defines 

f(t) = 2-e t, 
۱۱۱۳ ۱ = l and at the outset h(x,0) = g(x), and after one 
etre n(x,1) = (1.633)9(x). 


Assume the following conditional distribution for S(t) is 


given. 
So = EE ps (t) ۳ LE 6 
pat) = or p (t) = Rost 
pa (t) = _ہ‎ CE pj (t) = E 7 
pj4 (8) = e "۳-۵ °° pjg(t) 7 1-e^* 
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Figure 6. Coherent Systems Associated with the Absolute 


Effectiveness Function of Example 1. 
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By (3-8), 
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Dt) 168 a + Ep = ice e 
Br l  -6t 1 -7t 
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d in ERECTO 3 -5t 1 _-6t 
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fes chiese -3t 3 -4t 1 -5t 
ES) CO” = € + "e = ag e 
6 > - e ^* 7 TRE i = Ni 
E AE تن سل‎ 
I(8,t) = 2t - 3- E Das - = ni 
oy (3-7), 
8 
a(t) = ) ati) 7676 
i=l 
EESC E) GEI» 02 T SU 
+3 I(4 t) + 2 1 (5.,t) > rer 
+0 I(7,t) + 0 I(8,t) 
۳ 2 3 e 9t _ 6 art 7 ا‎ e Ot 
280 8 7 3 
ای جج ئا‎ MEE er: 
5 4 3 


Now the given interval probabilities are equivalent to 
the following reliability functions of the associated systems 


in Figure 6: 
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R4 (t) = e 
_ Dt 
R, (t) =O 
o -4t 
3 L 
Ry (t) = e 
o -2t 
R. (t) = e 
EE 3-10), 
Jut) = س‎ - => E + <> 7ئ‎ 
J(4,t) = = - + a a amis 
a -2t il -3t 
Jis) = 3 = E a oa 
And by (3-9), 


5 
q(t) (9541795) 1 
=] 


نا 


SU OR AE TF NO ۶۹ 7۸+ 


+ 1 J(4,t) + 1 J(5,t) 





Aa 2 = 
= 22 E تہ‎ e I + e A + = 
3 E l ie l -3t 


which agrees with the previously computed posce 
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IV. UNIT PERFORMANCE MODEL 


The quality function, q(t), has already been introduced 
de measure of effectiveness.. The next step is the introduc- 
tion of a model to translate the user's judgment and intuition 
about combat dynamics into a description of the stochastic 
performance process of the proposed unit structure. The pri- 
mary requirements for the model are that it be sensitive to 
Structural alternatives, and that it provide results in a 
form suitable for input to the quality function. There are 
two alternative forms of performance information which can be 
used in calculation of q(t). If expression (3-9) is used, the 
ENDHts required are the reliability functions of the associ- 
ated coherent systems. If the calculation is performed accord- 
mage tO expression (3-7), the 9 SE needed is the set of 
interval probability functions {p,j(t),....Pp}y(t)}. This 
Second type of information is the type of performance descrip- 
mon which is ost easily determined employing the model 
proposed in this section. 

The unit to be modeled is assumed to be a complex struc- 
ture which includes, as components, personnel, equipment, and 
ee mnels of information. The assumption will be made that 
each of the components is either funetioning or failed. “At 
the opening of an engagement all components are considered to 
be functioning. Some components such as personnel are non- 
renewable and do not return to the functioning condition once 


they leave it. Other components are renewable ana can 7 
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to the functioning condition repeatedly. For instance the 
information channel between a squad leader and a team leader 
which consists of voice communication can be interrupted by 
battle noise, then restored when the noise subsides.  There- 
fore an engagement can be looked upon as a series of events 
in time which change the state of the unit, through attrition, 
interruption, restoration or some combination of these effects. 
The state of the unit remains unaltered between these events. 
One approach commonly used to portray this process is to 
consider the performance of each component to be probabilis- 
tically independent. This assumption is not acceptable for 
this application since it severely limits the consideration 
of dependence in the attrition process. 

In order to integrate structural sensitivity into the 
model, it is necessary that the interevent times and the na- 
ture of each event be considered to obey random distributions 
which depend upon the surviving structure. The following 
Peeposal for accomplishing this is derived in part from the 
work of Marshall and Olkin [Ref. 9]. 

The engaged unit is subject to a hostile environment which 
consists of the delivery over time of one or more types of 
Heehal projectiles, including nuclear weapons, artillery, and 
small arms fire. Two assumptions are made about these fires. 

1) The delivery of each type of projectile constitutes 

an independent time homogeneous Poisson process, 
having interarrival times distributed exponentially 
with expected value - where k is the index of the 


k 
weapon type. 
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2) For each unit state and each weapon type the damage 
caused by each round obeys a probability mass func- 
tion: Pee) = P [Projectile of type k causes 

transition to state j | unit is 
in state i]. 
Lom these assumptions it follows that when the unit is in 
Ec 1, the distribution of arrivals of projectiles of type 
k, which send the unit into state j, is Poisson having expo- 


nential interarrival times with expected value PONTI 
k E 


Now in general the unit in state i can be subject to 
transition into state j as a result of several threats. If 
there are z types of projectiles which can cause the transi- 
tion, then the distribution of all such events is the super- 
position of z independent Poisson distributions, which is 
itself Poisson, having exponential interarrival times with 


1 


mean 
Py (9117, + ... + P_(3]1)t 


7 

In addition to attrition events there are also interruption 
and restoration events which describe alterations in condition 
of the renewable components. It is assumed that for each state 


i all functioning renewable components have an exponential time 


to failure with mean 3 where j is the state which the com- 
1] 


ponent's failure causes the unit to enter. Likewise, if state 


i has a failed renewable component, then = is the mean of the 
ik 


exponential time to restoration and consequent transition into 
State k. 
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Let as be defined as the time interval from the unit's 
entry into state i until an event occurs causing transition 


into state j. T is then the time until an attrition event, 


ij 
or an interruption-restoration event of the proper type 
occurs. This is again a superposition of independent Poisson 


processes resulting in a Poisson distribution having exponen- 


tial interarrival times with expected value, 


I | 1 


En E (۰ SR + lob ES + E x Aii. 
and survival function 
-(Ai.)t 
P[T.. > t] = e J ۲ 


1) 

In general the unit in state i faces threats which can 
cause transition to several different states. The outcome is 
determined by the type of the earliest state-changing event to 
occur. Let T, by defined as 'the length of the time interval 


between transition a-l and transition a, so that 


8 


0 
state i was entered at transition a-l. If transition from 
State i to state j is not possible, then De is considered to 
be infinite, with ipe = 0. Howard [Ref. 6] refers to this 

type of model as a "competitive process" since it is equivalent 
to a process in which, upon entry into state i, the random 


variables fia are sampled and, in effect, compete for the 


earliest outcome. 


41 


E B et mo To" سپ‎ PR a o vw 


rs 


a 





Rue survival function of T is derived as follows: 


۱۳ 5 ۷ ۳] Tea | 


0 2 72 ıN 
= PUT, 1 A] P[T,, > tle. -PITIN PNE] 
“Ag TE TE es Ain) t 


= e 
Thus T. is distributed exponentially with mean, 


il 
2 S EE LE PE 


e 


Thus the original assumptions of Poisson arrival of projectiles, 
and exponential interruption-restoration times, imply that the 


times between state-changing events are also exponentially 


distributed with mean = : 
Qi 

The validity of these assumptions must be considered in 
any application of the model. The assumption of Poisson inci- 
dence of fires is characteristic of relative stability in 
delivery rates. This distribution would not be consistent 
with scenarios which included large-scale enemy reinforcement 
or temporary massing of fires. 


The exponential distribution of the interevent times, T. 


implies the following result: 


Il 
Il 
(D 


Di ۴8 E 
a a 


PIT, > s] 


This result is referred to as the "memoryless" property or, in 


reliability theory, as the "no wear" property. The effect of 
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Me property in this model is that a unit's probability of 
successfully remaining in its current state for one addi- 
bwonal minute, is completely independent of the amount of 
time already spent in the state. This property fails to ac- 
count for such factors as fatigue or depletion of supplies. 
However in any case where the effect of these factors is 
insignificant compared with the actual attrition and interrup- 
tion, the "no wear" property may prove to be acceptable. 

The major consequence of the exponential assumptions is 


mut the Markov property holds? That is for t > لا‎ < 0, 


dc) = di) > v 290) = PIS po 


or the distribution of future S(t) depends only upon the cur- 
rent state, and not at all upon earlier information. The 
disadvantage of this property is that it does not permit con- 
sideration of such phenomena as momentum and collapse of 
situation. In applications where such effects are considered 
important, alteration would be required, 

The fortunate consequence of the Markov property is that 
it qualifies the process, {S(t), t > 0}, as a continuous para- 
meter Markov chain, a type of process Susie possesses a well 
developed method of solution. 

The decision-maker using this model is presumed to have in 
mind a specific unit structure and a specific engagement 
scenario. From the scenario he determines the variety of 
projectile types to which the unit e vulnerable along with 


mates of fire Tye Then for each state i he must specify the 
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mass functions EO) for each type threat k. These mass 
functions are the primary ingredient of the model's sensi- 
tivity to structure. They should be based on such variables 
as each individual's level of training, the capacity in which 
he is serving due to attrition of key personnel, the degree 
to which he must expose himself to hostile fire as a result 
of his current position, and the current condition of channels 
of communication. The specification of the mass functions 
will be the most difficult part of parameterizing the model. 
But it is indispensible if worthwhile investigation of unit 
structure is desired. 

The fact that the Markov property holds leads to the 
derivation of two sets of simultaneous first order differen- 
tial equations which can be solved for all the functions 
Pig tt). These two sets of equations are known as Kolmogorov's 
forward equations and backward equations. The actual deriva- 
tion is not presented here but is given in Ross [Ref. 11]. 
Before presenting these equations, several notational concepts 
are introduced. 

There are N^ distinct functions FE which appear in 


the functional matrix 


44 





Aa 2 ۲۳ he. o e eh ور‎ - mo 


where the function E appears in the ith row and jth 


Column. For any t 2 0 the rovs of D(t) add to one. Since 


p,:(0) = 1 is] 


J 
= 0 ifj, 
it follows that 
10 0 
ea ۲ 
P0) (= = ۰ 
; 1 0 
MES DI 


1۳ Bunet:ional matrix PUE) is defined as follows: 


p'ij(©) ... p'jy(t) 
P'(t) = 
Pee aE) او‎ E EE 
| Nl NN | 
where 
d p; (t) 
' = J 
و‎ de 


Using the parameters already defined in this section, the 


Sonstant matrix A is defined as follows: 


NI 
where E is the reciprocal of the mean of B and Aij appears 
ith row and jth column. Since no transitions are considered 


from any state into itself À ii A O NS 
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The constant matrix B is defined as follows: 


À Ax) He Ain 
Ay] 7 A) : 
B = $ ۰ ۰ 
NA Ay 1 
N 
where RE B. the reciprocal mean of the interevent time 
j=1 when the unit is in state i. 


Kolmogorov's equations are stated in matrix form below. 


The backward equations are given by 


(4-1) PO) = DEP OM 
with initial condition P(0) = I. The forward equations are 
given by 

(4-2) | PUE) ١+٤٤+٢) ٣ 
with initial condition P(0) = I. The operation indicated on 


the right of each expression is matrix multiplication. The 
fact that B P(t) = P(t)B is of special interest since matrix 
۲۰٠-۰ plication 1S not, in general, commutative. The solution 
of the process will consist of solving the backward equations 
using Laplace transforms. 

The Laplace transform of any functional matrix Alt), is 


indicated by 
۱ (C). 


which is simply the matrix of Laplace transforms displayed in 


the same order as in A(t). For any functional matrix A(t), 
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(2-3) L{A'(t)} = A'(s) = sA(s) - A(0), 


an equation proved by Kreyszig [Ref. 8]. Also for any func- 


wnal matrix A(t) and any constant matrix C, 
(4-4) Lie zo er 
(4-5) Tans) ME 


If Laplace transforms are taken of both sides of (4-1) 


the result is 
1:۱۳ ۹)9 ۳۲ 
which by (4-3) and (4-4) is equivalent to 
sP(t) - P(0) = B.P(t) , 
or since P(0) = I, 


sP(t) - B P(t) =I, 


OL 


~s 


isl = BI Pl SÊ, 


Premultiplying both sides by the inverse of [sI - B] gives, 


P(t) = fsı - Bf, 


which gives the solution for L{P(t)} in terms of the constant 
۳1 B. The solution is complete when the transform is 
inverted to give P(t). 

Howard [Ref. 6] shows that 


L(P(t)) = [st - B] 1, 


sand only if 


P(t) = ess f 
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where 


2 3 
م٣ خط + 1 ے‎ + Er +2 t 


is called the matrix exponential function. This expression 
may be used to calculate approximations to the elements of 
Hi), but will not provide exact solutions. Such a solution 
can be arrived at by finding the inverse of [sI - B] and 
using the method of partial fractions to translate the result 
into invertible Laplace transforms. 

Example 2 

Consider again the tank killer team of Example 1, with the 
same three components, all considered nonrenewable, and the 
me state function Six). 

The scenario of interest involves two types of external 
en ac. Tne finst threat consists of Small arms pro; ceerics, 
delivered with a Poisson rate of fire o 50/hour. The 
second threat is artillery fire with rate a VO nour: 

The decision-maker now determines the approximate proba- 
bility mass function for each state-threat combination. 

For state one and the small arms threat, the user first 
notes a single round cannot cause more than one casualty. He 
maerore defines Be) = ۳ E P, (7/1) E P, (8|1) R 
since these transitions would entail multiple attrition. He 
feels that the gunner is about 50% more vulnerable to small 
arms than the leader because of his firing position, and the 
Fact that the backblast of his weapon pinpoints his location. 
The weapon is considered only half as vulnerable as the leader 
because its small presented profile. The user thus estimates 


the following mass function: 
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POD P (5|1) = 0 
۳۰ ۱۶۳۳ . 5 Pec 
P, (3/1) = .02 P,(7|1) = 0 
ئ۶‎ + 7 701 77 “ “) ۱۶۶۶۵٣ 


where the actual magnitudes may be based on field data, 
Emulation, or intuition. 

For state one and the artillery threat multiple attrition 
is possible. The user feels that each artillery round has a 
capability to destroy all three components and defines 
p, (8/1) x Ec The gunner is considered more vulnerable to 
artillery fire for identical reasons as cited under small arms. 


The decision-maker defines this mass function as follows: 


P_(1|1) = ٤ ES ٤ 
70 2922.05 P,(6|1) = 0 
pone Po(7|1) = 0 
Pp, (4|1) = 5 Pr 00 MN 


The probability of transition into state seven is considered 
to be zero, since a round placed well enough to eliminate both 
personnel would be expected to destroy the weapon as well. 
Transition to state six is similarly considered impossible. 
The actual magnitudes are again based on the user's judgment. 
Since there are no renewable components the parameters = 


can now be computed: 
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L3 


14 


ES 


16 


7 


18 


Il 


0 (trivial transitions are ignored) 
P, (2|1)ty + P_(2|1)t, 


OSO 


+ 


(205) (20) 


NO مر‎ S 


۱۰ ما‎ ۱ Ca oo) 


paucae. + R 


(01) (50) FE COSCO) 


BCE A Te) ۱ li. 
(0) SOHC COO) 

1 

۶8 ۶۲ 9 CC ula. 
0 

p: BONUM 
0 

P,(8[1)1, + Bars 


0 + (.1) (10) 
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The construction of the full set of mass functions 
results in the constant matrix A shown in Figure 7, and the 
functional matrix [SI - B] in Figure 8. 

The method of solution is to invert [sI - B] by the method 


of cofactors using the relation 


(-1) 273 


PSP Re ۳. 
where |sI - B| is the determinant of the matrix [sI - B], and 


Mi; (S) is the determinant derived from |sI - B| by deleting 


the jth row and ith column. 
The advantage of the cofactor method is that it permits 
the elements of [SI - ES to be computed individually saving 


time in cases where only a portion of the inverse is needed. 


The results are as follows: 





= 2 -7t 
Pi, (8) = = Le |. 

p.n (s) = - = "۵٣ 
P12 S+7)(S+5 s+5 s+7 
= Lfe >t AA 
p.~ (s) = 2 MEM e A 
P13 (s+7) (s+6 Sd S47 

ee n 
۳ (s) = 1 = 1l = E 
P14 — (s*7) (s*6 8 +6 S47 
- Den 0 


Du 





0 0 0 2 2 
0 0 0 l 0 
0 0 0 0 ili 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
Figure 7. Constant Matrix A 
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S45 0 
0 s+6 
0 0 
0 0 
0 0 
0 0 
0 0 
Figure 8. 


zi -1 0 0 - 1. 
0 -2 -2 0 -1 
0 zu 0 ES 2 

510 0 -1 23 7) 
0 s+4 0 0 — À 

| 
0 0 s+4 0 EA 
0 0 0 s+4 -4 
0 0 0 0 7 
E 
Functional Macrix by 
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2 = 1 4 2 
و 1ظ‎ ( ^ (S47) (st+4 5 S+7) (s+5) (s+4 N s+7) (s+6) (s+ 
2. cas 2 
S+7 s+6 s+5 s+4 
Sen T Deane 
~ ۲ 4 
A EI A (stoned 
l 3 
= E = 2 - 2 + 2. 
S+7 S+6 S45 s+4 
= em 5 یط‎ >-6 7 Ja >t + E 
2 2 
p (s) = 2 
H7 (s+7) (s+6) (s+4 
2 3 
= EN = 2 + 2 
S+7 S+6 s+4 
T É e Dar ES 
2 2 
> Cu 48 ۳ 32 ۳ 6 
P18 E s+7) (s+6) (s+4)s S+7) (5+5) (5+4)5S s+7) (s+6)s 


+ 2 + 4 + - 1 
(s+7) (s+5)s (s+7) (s+4)s (s+7)s 


Lei O 2 
dou S+7 S46 S45 s+4 
= nen ae nie Oa ese ae 989۲٦ 
8 
As required 3 ES = +,4 OUR 
jj 


The interval probabilities {py Wr... Pig (t)) are now 
known and the quality function can be computed for any 


Functions f(t), and g(x) or dli).|| 
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It should be noted that the example was particularly 
simplified by the fact that the matrix [sI - B] had no non- 
zero entries below the diagonal. The reason for this is 
that all three components were considered to be nonrenewable. 
Thus when the unit left any state it could never return to 
it. And by arranging the states in order so that the earlier 
states corresponded to the higher rows, the triangular form 
was assured. The resulting denominator of all the expres- 


sions P. (9) was Of the form 


8 
|sI - B| = 0 ۱ 


i=l 


which is already in the factored form necessary for the partial 
fraction expansion. 

If however at least one of the components had been renew- 
able, there would: have been at least one positive entry below 


the diagonal and the denominator of the expressions p 


1318) 
would have been of the form, 
8 
isi = Bl = | ea) re 
i=l | 


where K(s) is a polynomial function of s of degree less than 
eight. Therefore |sI - B| is a polynomial of degree eight 
whose roots and factors are unknown in general. This computa- 
tional problem is formidable since the degree of polynomials 
to be solved will equal N = Di (where n is number of system 
components), and there is no exact solution for polynomial 


equations of degree greater tren folder. Mere 6 ٣ 


DI 








computer routines which perform approximation procedures, 
which may prove to be of value in employing this model. 

Finally, a rough measure of effectiveness may be computed 
from q(t) without inverting the Laplace transforms when 


mee) = 1, € 20. Let the quality Inde. Be EINE ات‎ ORONS. 


Q = q(%) = Las Si P;; (u) du 


1 


N ~ 
), 8G) By 5 (0) 


Thus Q is the expected aggregate value over an infinite time 
period. Of course, any realistically defined unit will 

reach a state with zero effectiveness in finite time, so 

that Q has a finite value. However this measure does not dif- 
ferentiate between units whose quality functions are thinly 
distributed over the real time line and those whose effec- 
tiveness is heavily concentrated in short intervals of reason- 


able engagement duration. 
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Vo CONCLUSION 


The purpose of this thesis has been to outline a proposed 
modeling viewpoint. It has not been advanced as a finished 
product. It is hoped that the concepts introduced will be 
considered as a worthwhile point of departure for additional 
research, which could result in establishment of a valuable 
Ronson Tool for choosing among alternita n TIC ELLES 

The most critical unsolved problem is the computational 
difficulty involved in solving the model. In addition to the 
inversion of Laplace transforms, the current solution neces- 
sitates the evaluation of large order determinants, which can 
become unwieldy even with electronic computers. Thus the 
current span of application is limited to small organizations. 
Research to extend this span should concentrate on finding 
Simplified computational procedures for exact solution, or 
alternatively, determining efficient and accurate approximation 
methods. The use of the infinite series definition of the 
matrix exponential function as an approximation method was 
mutroduced in Section IV. Unfortunately, this expression 
requires sequential multiplication of large order matrices, an 
operation which is as computationally explosive as evaluation of 
determinants. 

Current research in reliability theory may offer a valuable 
approximation method, particularly since ali of the random 


processes and measures of effectiveness were shown to be 
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expressible as linear functions of variables related to 
reliability of coherent systems.  Esary and Proschan [Ref. 4] 
show the usefulness of reliability bounds and approximations 
based on the concepts of minimal cut sets and minimal path 

sets of components. Adaptation of this method to the avail- 
ble type of data inputs may prove useful to Ene current modela 

A possible exact solution method using flow graph analysis 
is proposed by Howard [Ref. 6], but at present it appears to be 
subject to an even greater computational burden than the cur- 
rent method. 

In order to enhance the accuracy of the model it may be 
desirable to alter the assumption of exponentially distributed 
interevent times in favor of some other distribution. This 
meula change the process to one of me sen Markov va E E 
These processes are more difficult to solve in general, but 
offer greater ee in definition. Ross ١ 57 
Howard [Ref. 6] consider several solution methods. 

A highly desirable addition to the model would be a provi- 
sion for random threat. Rather than assuming constant Poisson 
rates of enemy fire, such a refinement would allow rates to 
vary randomly, thus representing periods of non-engagement as 
well as change of situation during engagements. The analogous 
reliability concepts of random wear and random shock have been 
studied by Gaver [Ref. 5] and Reynolds and Savage [Ref. 10]. 

Despite its capacity for refinement, the proposed model in 
present form does fill many of the needs for the study of struc- 


Lure. It offers the decision-maker a model which uses his own 
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judgment, in an intuitively clear process of parameterization, 
in order to compute a functional measure of effectiveness 
based on fundamental functions which he, himself, defines. 
Indeed, with the exception of clearly stated assumptions made 
in building the model, and the purely mathematical solution 
process, the model in use is completely established by the 


user, with whom ultimate decision responsibility rests. 
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